Poincare recurrence and measure of hyperbolic and nonhyperbolic chaotic attractors 
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We study Poincare recurrence of chaotic attractors for regions of finite size. Contrary to the 
standard case, where the size of the recurrent regions tends to zero, the measure is not supported 
anymore solely by unstable periodic orbits inside it, but also by other special recurrent trajectories, 
located outside that region. The presence of the latter leads to a deviation of the distribution of the 
Poincare first return times from a Poissonian. Consequently, by taken into account the contribution 
of these recurrent trajectories, a corrected estimate of the measure can be provided. This has wide 
experimental implications, as in the laboratory all returns can exclusively be observed for regions 
of finite size only. 



In dynamical systems one is often interested in calcu- 
lating asymptotic invariant physical quantities that are 
independent of the initial conditions and invariant under 
the evolution of the dynamics. Those invariant quantities 
are the main subject of ergodic theory, which was suit- 
ably applied to chaotic dynamical system, producing im- 
portant average estimates, like Lyapunov exponents and 
dimensions jlj]. Chaotic systems are ergodic and there- 
fore time averaged quantities calculated from a typical 
trajectory can also be calculated through space integrals 
using the natural measure of a chaotic attractor, which is 
related to the probability of finding a typical trajectory 
in some region of the attractor. Hence, the calculation 
of this measure yields essential observable quantities of 
chaotic systems. 

It is well accepted that the support of the measure as 
well as properties of the attracting set are hierarchically 
approached by the set of unstable periodic orbits (UPOs) 
which are embedded in the chaotic attractor j^i^ JJ. For 
hyperbolic systems this approach allows a reconstruction 
of the fractal dimension of a chaotic attractor through 
UPOs 12, For the important class of nonhyperbolic 
chaotic attractors, i. e., when there exist UPOs whose 
stable and unstable manifolds exhibit tangencies, only 
numerical results for averages in phase space exist, sug- 
gesting that, on average, the nonhyperbolic regions do 
not play a special role ^ . 

Recently, an exact result on the probability distribu- 
tion of the series of first return times (FRTs) Ti,i — 
1,...,N for measurable dynamical systems (including 
nonhyperbolic ones) was provided . It was proven that 
the distribution of the FRTs, if the return time is larger 
than some constant value and the size of the region goes 
to zero, approaches the Poissonian p{t,B) = /zexp"^'^, 
with /J, being the probability density measure of the re- 
gion B of size e, as also discussed in [^Q. Thus, 

limp(T>^,S)-Aiexp-^-<G(M^)) (1) 

with C being a suitable normalizing factor and G{fi{B)) 
a function of the probability density measure, fx, in B. 
The multifractal spectrum of return times has been in- 
vestigated as well [lOl |. 



In this work, we extend this rigorous treatment to the 
case for which e is far from zero and r is arbitrary. This 
is motivated by the fact that experimentally only returns 
to regions of finite size can be observed. In the following 
we take the returnung region to be a square of size e in- 
tersecting the invariant set (attractor) . We show that for 
both, nonhyperbolic (logistic map and Henon map) and 
hyperbolic (Cat map) systems, the presence of recurrent 
trajectories that are not associated with UPOs inside the 
region increases as one increases the size of the region, 
causing the measure to become abnormally singular. In 
the case of nonhyperbolic systems, even for square of very 
small but finite size, the contribution to the measure of 
these recurrent trajectories is already large if the interval 
is close to a homoclinic tangency. Hence, as one observes 
a dynamical system for intervals of finite size, the local 
dynamics in B can neither be completely governed by the 
linearization of these UPOs as proposed in the theorem of 
Hartman-Grobman llj, nor the measure inside the inter- 
val, fJ.{B), can be exactly calculated through the UPOs 
inside it. However, from the distribution of the FRTs 
we show that it is possible to calculate the measure ex- 
actly by introducing a correction of the measure from the 
eigenvalues of the UPOs inside the interval, when the in- 
tervals are far from small. 

As already mentioned, an important result concerns 
the reconstruction of the measure using UPOs 5j. For 
hyperbolic systems, one can derive a formula to calcu- 
late the probability density, fi{B) eig of small squared 
subregions B in the attractor. More explicitely 

= (2) 

where the L^. are the positive eigenvalues of all fixed 
points located in B of the j-fold iterate of the map, (i. 
e., the fixed points are period-j UPOs G B), and Nj rep- 
resents the number of period-j UPOs e B. 

We conveniently define a finite sized region to have a 
hyperbolic character if that region exhibits a distribution 
of FRTs close to a Poissonian. In this case, the measure, 
as well as the dynamics in that region, is described by 
the UPOs inside it, and typically there lie neither homo- 
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clinic tangencies nor low-period UPOs inside the region. 
Analogously, we define that a region is nonhyperbolic 
if the distribution of the FRTs deviates from a Poisson 
law. In that case, recurrent trajectories play an impor- 
tant role in the dynamics and in the contribution to the 
measure; here, typically low-period UPOs or homoclinic 
tangencies are found. We show that the larger the con- 
tribution of the recurrent trajectories to the measure is, 
the larger is the deviation of the calculation of the mea- 
sure from the eigenvalues of the UPOs inside these in- 
tervals. The contribution becomes negligible as either 
the interval size tends to zero or the intervals are placed 
in hyperbolic regions and off homoclinic tangencies and 
low-period UPOs. 

The probability measure in a square B due to all the 
recurrent orbits (including the UPOs) with recurrent pe- 
riods between Tfc and t; can be calculated from the dis- 
tribution of the FRTs p{Ti, B) as 



(3) 



where the average first return time is given by 

1 " 

{r{B))= lim - Y.n{B). (4) 



Note that if one sets in Eq. ||2Jl tj, = Tmin (minimal 
FRT in B) and t; = Tmax (maximal FRT in B), then, 
X^ii" p{Ti,B) = 1 and one recoveres Kac's theorem, 
which states that 



/i(i3) 
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with {t{B)) the average FRT, relating a time averaged 
quantity {t(B)) with a spatial quantity, /i(B). 

In order to calculate the measure exclusively due to 
UPOs inside B, we hypothezise that the distribution of 
the FRTs can be split into two discrete functions (for 
a small FRTs) and into two continuous functions (for 
large FRTs), one describing the contribution due to the 
UPOs and the other the contribution of the recurrent 
trajectories, respectively. Hence, 



M'(^) = m%.EC + mUPO + mREC + KBTUPO, 



with fi{B)%Ec = (7^^»=-«Ec^(^''^)' ^(^);7Po 

J^J2^=ruPo and 
1 



(6) 
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/?fl£cexp(-"^)dT (7) 
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{r{B)) 



a exp 



dr. (8) 



The index d designates a discrete summation and c in- 
dicates a continuous integral. fJ'{B)j^^(j is the measure 



due to short recurrent trajectories, fi{B)^pQ due to low- 
period UPOs, fiiBYi^EC measure due to long recur- 
rent trajectories and ijl{B)'ijpq due to high-period UPOs 
'l4'|. The return time trec stands for first returns that 
are different from any of the period of the UPOs inside B, 
and Tupo for all periods of the UPOs. The time r™pQ 
corresponds to the minimum period for which one can 
consider the distribution of FRT to be well described by 
a continuous exponential function and the measure can 
be calculated with the Eqs. and 0. These two in- 
tegrals are constructed under the assumption that the 
probability distribution is given by 



p(t, B) = (3 exp " 



(9) 



i.e., it agrees with Eq. iQ. Since ^{B) cx e^", A, > 
being the pointwise dimension for that interval 13], as 
e — > and t > ^ , the distribution © should approach a 
Poisson function p{t,B) = /icxp^'''^. Equation ^ can 
be broken in a summation of two distributions, one due 
to the long recurrent trajectories p'j^EC ~ Prec exp(~"'^\ 
and the other due to the large-period UPOs, a Poissonian 
of the type Pijpo = aexp^""'^-', such that p'pEC + Pupo 
= (3exp^~°'^\ Therefore, /3_rec + a — P, and, from Eq. 
(0), as e — > 0, Prec — 0, and a — p,. The coefficients 
f3 and a are obtained by fitting the distribution of the 
FRTs by an exponential of the form of Eq. 1^. Inte- 
grating Eqs. O and (jSJ we get Pp^c (a ~ "fP 



exp" 



Pupo = with 7 = exp ' 

Assuming that Pp^c ^^'^ Pupo ^^^^ negligible in com- 
parison with PpEC a^*^ Pupo^ thai p^eig = Pupo 
(what is true for moderately nonhyperbolic regions which 
are predominant in the attractor), we arrive at a formula 
that expresses the value of the measure through the mea- 
sure calculated exclusively with the contribution from the 
UPOs inside B 



fi{B) - Kfi{B)EiG, 



(10) 



where 



Our model of a nonhyperbolic system is the Henon 
map, H : X ^ X, x^+i = a — + byi, and tji+i = Xi, 
with a = 1.4 and h — 0.3. The nonhyperbolicity of this 
map is due to tangencies of the stable and unstable 
manifolds of periodic orbits embedded in the chaotic 
attractor. 

In Fig. n] we show what should be expected from a 
hyperbolic and from a nonhyperbolic region. In Fig. 
n^a), points represent all the UPOs up to period 23, 
calculated using the method of Ref. j^], and stars 
depict the primary tangencies together with their 
five images and preimages, where primary indicates 
that the curvature of the manifold is minimal in their 
neighborhood. In Fig. Q^b) we display a typical nonhy- 
perbolic region (box) centered at the primary tangency 
Ti : {x,y) = (1.7801,-0.0949) with box size e = 0.02. 
The lowest periodic orbit found in that region has period 
18, whereas the smallest FRT is r = 9. Points represent 
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the Henon attractor, and the filled circle a component 
of a period-9 orbit. The stripe marked by Ws pictm'es 
points inside B that under (the 9-fold iterate of the 
map H) remain in the interval, represented by Wu- As 
one can see, the stripes Wg and Wu are aligned along 
the stable and unstable manifold of the period-9 orbit 
located outside B. Consequently, the tangency creates 
recurrent trajectories that are not associated with any 
UPO inside B. Note that, as could be anticipated for a 
nonhyperbolic region, the stripes Wg and are almost 
parallel close to Ti. In Fig. QJc) we demonstrate that 
a similar effect occurs when there is a low-period UPO 
inside B (here of period-2). Although there exists no 
period-13 UPO, we nevertheless find a FRT of r=13. 
That is caused by points on the stripes Wg and Wu, 
which lie on the stable and unstable manifold of a 
period-13 orbit outside B. Finally, we show in Fig. ^d) 
a typical hyperbolic region, where to all observed FRTs 
of length T a period-r UPO inside B corresponds, and 
the manifolds cross transversally. 



returns to that interval as before. This gives with Eq. 
(jHJ /i' = 7.2295 X 10~*, whereas the exact measure, 
using Eq. ensues /i = 7.2422 x 10""*, very close to 
the value obtained with the correction formula Eq. 10. 
Again, most of the measure is due to the UPOs. From 
Eq. the best estimate (calculated for UPOs of up to 
period 30) of (Ieig = 7.5312 x 10"^. For this case, Eq. 
(|1(J|I has no use since heig > 
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FIG. 1: Points represent the chaotic attractor F, filled circles 
periodic orbits, and Wa denotes points inside B that under 
(where P is the period of the UPO) remain in the inter- 
val; their P-fold iteration is shown by Wu- In (b) the box 
is centered in a primary tangency and in (c) in a low-period 
UPO. Both are typical nonhyperbolic regions. In (d) a char- 
acteristic hyperbolic region is depicted. 

To quantify now these findings with our theory, we 
analyze two specific intervals with e — 0.02, a nonhy- 
perbolic one, Bi, centered at the primary tangency Ti, 
see Fig. Qlb), and a hyperbolic one, B2, centered at the 
point {x,y) = (1.1181,0.1472), see Fig. Hid), for which 
the smallest period of all UPOs is found to be 21, and 
no tangency is present. 

For the nonhyperbolic interval Bi, we get 
H%EC — 1-28 X 10"'* from the return times trec = 
(9,16,17,18,19,20,22), <po=0, fi^^.^ = 2-41 x 10-^, 
and ^l1Jpo = 2.885 x 10-^, with r^J,-^ = 23 and 
Tmax — 4284 in Eqs. JHJ, obtained by measuring 

400,000 returns to that interval. Employing Eq. 10, 
this yields a measure of /i' = 3.0923 x 10""^. We stress 
that the measure evaluated using Kac's lemma with 
400,000 returns, considerered to be the exact one, leads 
to /i = 3.100 x 10^"^, very close to our estimate from 
Eq. ©. We see further that most of the measure is 
due to the UPOs and the contribution from rccurcnt 
trajectories fijiEC I^rec neglegted. From 

Eq. (|2Jl the best estimate (determined for UPOs of 
up to period 30) is fiEiG = 3.003 x 10"^. Note that 
fiEiG — f^BOTH + l^upo- spite of this region being 
strongly nonhyperbolic, we apply nevertheless the 
correction in Eq. ((TUI), which yields /i(S) = 3.027 x 10"^, 
a value closer to the exact ^ than /j, eig ■ 

The hyperbolic region B2 results in nfiEC ~ ^' t^upo ^ 
0, = 9.97 X 10-^ and fi'fjpo = 7.712 x 10"*, 

with T™" = 23 and r^ax = 16871 obtained for 400,000 



FIG. 2: The weighted measure tiEllRM^^ i=i, 2, due to recur- 
rent trajectorie with respect to e for the interval Bi (circles) 
and B2 (squares). 

To understand how the recurrent trajectories 
contribute to the measure of these two intervals 
when the boxsize e varies, in Fig. [5] the value of 
^^^'^^^'^ = (a - l)Ai(Si),i = 1,2 is plotted against e. 
It can be calculated using only the information of the 
FRTs, and thus there is no need of knowing t^^Iq , a nu- 
merically envolving task. We see that for the interval Bi , 
centered at the primary tangency Ti , the contribution of 
the recurrent trajectories HREciBi) is much larger than 
i^REc{B2) for the interval B2, provided e > 0.0002. It 
grows moreover much faster with increasing e. Also, for 
B2 this contribution decays smoothly as one decreases 
e. This smooth decay in arbitrary intervals of finite 
size takes typically place in hyperbolic regions when 
the size appoaches zero. It is also encountered in the 
hyperbolic cat map as well as in the hyperbolic regions 
of the logistic map 

In Fig. m we plot |((4|^))| (squares), 
\{{Bm^J^mMlo)\ (circles), and |(( V(«.)- W.)..o ^| 

(diamonds) versus e, in a log-log graph. The average 
is performed over intervals Bi centered in consecutive 
points of a trajectory of length 5000 and fJ,{Bi)EiG is 
calculated using the set of all UPOs of period 29. The 
standard deviation bar is also shown in this figure. 
Whenever ^ < 0, which indicates the presence of 
low-period UPO inside the interval, fJ-{Bi)fjpQ connot 
be neglegted, and therefore the approximation proposed 
in Eq. (|10|l cannot be used. Furthermore, the calcu- 
lation of the measure from Eq. (O oscillates widely 
in intervals dominated by low-period UPOs, since one 
considers different sets of UPOs with varying periods. 
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FIG. 3: Averages (squares), 

|^( .(B.)-M(B.)..o )| (circles), and |(( m(B,)-a-m(B.)^xo ) | 
(diamonds), i = 1, 2, 5000, against e, in a log-log graph. 

Consequently, the averages calculated in this figure are 
restricted to intervals with ^ > 0. Additionally, r™pQ 
is considered to be approximately given by t^ec- 
a result, in an average sense, the presence of recurrent 



trajectories (circles) significantly affects the correctness 
of the calculation of the measure from the eingenvalues 
of the UPOs in Eq. ^ (squares). It is apparent from 
the figure that, on average, the correction proposed 
in Eq. for ^jleig yields a value closer to the real 

measure /i, as well as a lower standard deviation. 

In conclusion, finite size regions of chaotic attractors 
can generally be classified in two categories, hyperbolic 
and nonhyperbolic ones. For the former, the measure 
is almost completely supported by the UPOs inside the 
region and the distribution of the returns close to a 
Poissonian, while for the latter, the contribution of the 
recurrent trajectories not associated to UPOs inside the 
region is significant and the return time distribution 
deviates from a Poissonian. In this case an exact 
calculation of the measure can only be performed from 
the first return times. This has strong implications on 
experiments, as there only finite region can be monitored 
and many systems are nonhyperbolic. 
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